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It is demonstrated using three-dimensional computer simulations that some simple non- 
interacting electron models that include electron scattering by grain boundaries, exhibit 
coexistence of large persistent currents and small conductances, similar to that observed 
experimentally in isolated micron-scale gold rings. Models with simple grain boundaries, and 
models with small numbers of regularly stepped or atomically rough dilute grain boundaries have 
been studied and found to yield similar results, which differ markedly, however, from the 
predictions of models that assume only random impurity scattering. This difference is due to the 
fact that equilibrium persistent currents and non-equilibrium transport coefficients are physically 
different things and depend in different ways on the topology of the defect structure in a conducting 
ring. Experiments on metal and semiconductor rings that should be able to determine whether this 
is the explanation of the effects observed by Chandrasekhar et al. (Phys. Rev. Lett. 67, 3578 
(1991)) are proposed. 
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1. Introduction 

When a conductor in its ground state is immersed in a magnetic field, it acquires a magnetic 
moment and an associated circulating electric current. This current, being an equilibrium property, 
does not dissipate. It is thus referred to as a "persistent current." The basic theory of persistent cur- 
rents in normal metal rings was formulated by London^ in the 1930' s in the context of aromatic 
molecules. It was developed further by several authors in the 1960's.^~^ More recently, Biittiker, 
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Landauer and Imry showed that these currents should occur even in the presence of disorder. Che- 
ung, Riedel and Gefen^ then estimated the magnitude of the persistent currents that should flow in 
small, disordered normal metal rings with non-interacting electrons. They predicted that, at low 
temperatures, the typical persistent current should be I^yp ~ Iq II L where Iq = evpIL, e is the electron 
charge, L is the ring's circumference, is the electron Fermi velocity, and / is the elastic mean 
free path. Since then, persistent currents have been observed in an array of 10^ micron-scale copper 
rings by Levy et aZ, and in individual gold and semiconductor rings by Chandrasekhar et al 
and Mailly et al. The measurements on the single gold rings presented a puzzle that has re- 
mained unresolved: The measured persistent currents were of order /g in rings for which II L -0.01, 
i.e., the observed currents were two orders of magnitude larger ihaa had been predicted^ for rings 
with such small transport mean free paths and correspondingly small conductances. 

During the last three years, there have been several attempts to explain this observation as a 
novel effect of electron-electron interactions: Very large enhancements of the persistent current 
due to many-body effects were predicted in three-dimensional rings^^'^'* and also in one and two- 
dimensional models. However, subsequent more detailed three-dimensional calculations ' ' 
yielded no significant many-body enhancement of the persistent current. Also exact results of Lut- 
tinger liquid theory for defect-free rings as well as exact numerical solutions of models of small 
systems with and without defects^^"^^ have indicated that electron-electron interactions do not 
yield any major enhancement of the persistent current in one-dimensional rings. A recent theoret- 
ical study of two interacting electrons in a narrow two-dimensional ring also showed no enhance- 
ment of the persistent current by electron-electron interactions. Furthermore, the results of 
persistent current measurements on two-dimensional, few-channel, quasi-ballistic semiconductor 
rings^^ have been in reasonably good agreement with the predictions of the non-interacting elec- 
tron theory,^ so that it seems that there is no large many-body enhancement of the persistent current 
in those systems. It is also important to realize that every geometrical dimension of the gold rings 
in the experiments was very much larger than the electronic screening length, and that the plasma 
frequency of gold greatly exceeds the characteristic energy scales that appear in the physics of per- 
sistent currents. Thus from the point of view of standard many-body theory, one would expect the 
gold rings to behave as a classic normal metal, and there seems to be no clear physical basis for a 
giant many-body enhancement of the persistent current. In view of the above, although the search 
for a many-body explanation of the experiments of Chandrasekhar et al^^ is far from over,^^ it 
seems reasonable to begin considering other possibilities. 

The purpose of this article is to demonstrate that some simple three-dimensional models ex- 
hibit the combination of large persistent currents and small conductances that has been observed 
experimentally in the gold rings. These models do not introduce any unconventional electron-elec- 
tron interaction effects in order to obtain large persistent currents. They are based instead on well- 
known materials and transport properties of thin gold films. 

The previous theoretical work has been on models in which random impurity scattering is 
the mechanism determining the electron elastic mean free path /, and limiting the size of the per- 
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sistent current. However, the rings of Chandrasekhar etal were fabricated from gold films grown 
on oxidized silicon, an amorphous substrate. Such films are polycrystalline, and experimental stud- 
ies have shown^^'^^ that grain boundary scattering dominates the electronic transport properties of 
thin polycrystalline gold films at low temperatures. Random impurity scattering and surface rough- 
ness scattering are less important. It will be shown below that certain simple models that include 
grain boundary scattering exhibit the combination of large persistent currents and small conduc- 
tances that is observed experimentally, but is not found in the random defect models. The large dif- 
ference between the results obtained from the grain boundary and random impurity models is due 
to the fact that the persistent current and conductance of a ring are physically different quantities 
and depend in different ways on the topology of the defect structure in the ring. 

It should be stressed that the puzzle presented by the experiments of Chandrasekhar et al^^ 
involves not only the persistent current but also the conductance of the ring, and that both of these 
quantities must be calculated on the same footing if the experimental results are to be explained 
convincingly. However, most of the previous theoretical studies have been devoted to the question 
whether a many-body enhancement of the persistent current is possible; much less attention has 
been given to calculating the conductance. In the present work both the persistent current and the 
conductance are calculated numerically, and both calculations represent exact solutions of the 
models considered. The reason why the grain boundary scattering models with non-interacting 
electrons that are discussed in this work exhibit large persistent currents together with small con- 
ductances, while random defect models do not, is due as much to the different behavior of the con- 
ductance as to the different behavior of the persistent current in the two classes of models. 

Grain boundary scattering is the dominant scattering mechanism at low temperatures in thin 
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films ofgold,^''^°and some other metals. It therefore needs to be considered in any quantitative 
explanation of the experiments of Chandrasekhar et al. ^ ^ Although, as is shown in this article, grain 
boundary scattering offers a possible explanation of the experiments of Chandrasekhar et al, 
whether this will in the end turn out to be the whole story is an open question. The reason is that 
the gold rings in question have grain boundaries that are widely separated on the atomic scale and 
this "diluteness" of the grain boundaries is very important. Three-dimensional metal rings with 
large numbers of widely spaced grain boundaries are difficult to simulate numerically because of 
the practical limitations of computers. Ab initio calculations of the electronic structure of such sys- 
tems are at present out of the question. Thus the present study has been limited to very simple mod- 
els with separable Hamiltonians describing rings with arbitrary numbers of grain boundaries, and 
to more realistic (but still simple) non-separable model Hamiltonians describing rings with small 
numbers of dilute grain boundaries. Some important unresolved issues are how many grain bound- 
aries the rings of Chandrasekhar etal^^ actually contained, and if this number was large, whether 
the present results can be applied to that case. It seems probable that these and many other ques- 
tions will need to be resolved experimentally, and several good experimental tests of the present 
explanation of the observations of Chandrasekhar et al ^ ^ are suggested at the end of this paper. 

It should be noted that the present paper addresses the problem of the typical persistent cur- 
rent in a single metal ring, the quantity measured experimentally by Chandrasekhar et al}^ The 
problem of the average ring persistent current over a large array of rings that was measured by 
Levy et al}^ although closely related, is significantly different. This is because the persistent cur- 
rents carried by different rings in the array can differ in sign, which results in strong cancellations 
in the average. The theories of the average ring persistent current in large arrays of rings that have 
been developed to date^^' have assumed that the electrons are scattered only by random impuri- 
ties. In view of the results for single rings that are presented in this article, theoretical studies of 
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arrays of rings based on models with grain boundaries would clearly be of interest, but are beyond 
the scope of the present paper. 

2. Grain Boundary Models and Method of Solution 

Consider a thin, multi-channel ring threaded by a magnetic flux, described by a 
Hamiltonian 

^ = - ^j,k,l + 

where C, = exp(27r/^>/Z), ^> is the magnetic flux threading the ring in units of the flux quantum 
h/e. Zis the circumference of the ring in units of the lattice parameter, at^^ is the electron creation 
operator at site jkl; the spin index is suppressed. The first term in // represents electron hopping in 
the azimuthal direction around the ring; the other two terms describe hopping in the two orthogonal 
directions (see left inset, Fig.l, for the hopping directions a, P and y). For a defect-free ring, the 
hopping coefficients tjj^i in H are all taken to be equal, tjj^i = tQ. Grain boundaries are modelled by 
setting those tji^i that represent hopping across a grain boundary equal to with tgi, < tQ. 

This is clearly a greatly simplified model of a metal ring with grain boundaries. However, it 
allows one to examine for the first time the effect on the persistent current of extended defects that 
cross the ring and partition it into "grains," as opposed to random impurities which do not. This 
topological difference between random impurity and grain boundary models is crucial, as will be 
shown below. Another important advantage of this model is that it can be solved numerically for 
both the persistent current and the conductance without making approximations. Thus all of the nu- 
merical results presented in this paper constitute exact solutions of the model. 

The persistent current / was evaluated at zero temperature by finding the eigenvalues of H 
numerically, and using the result^ 

I=-elhdEld^, 

where E is the total electronic ground state energy of the ring. 

In addition to the persistent current, it is essential to calculate the resistance R of the ring, 
since the value of the key experimental parameter IIL was determined from resistance measure- 
ments. For definiteness, let R be the resistance (in the absence of magnetic fields) of the wire made 
by severing the ring; the length of the wire is equal to the circumference L of the ring. In this work 
the resistance R is evaluated from the Landauer conductance formula 

G = le^Tlh, 

where R=l/G and T=Tr(tt^) is the total multi-channel transmission probability through the wire at 
the Fermi energy.^ ^"^"^ On the other hand, in terms of the effective mean free path, the conductance 
is 

G = S^e^l/{3hL), 

where r\ is the number of conducting channels at the Fermi energy. Comparing these expressions 
for the conductance yields l/L - 3T/4r[. Thus, in the present work. 
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will represent the quantity l/L which was inferred experimentally from measured resistances. 
The parameter T/r\ that appears in X has a simple physical meaning: it is the probability that an elec- 
tron at the Fermi energy is transmitted right through the wire made by severing the ring. In this 
formulation, the puzzle raised by the experiments of Chandrasekhar et al^^ and to be addressed 
here is why gold rings with low electron transmission probabilities (values of ?i ~ 0.01) are ob- 
served to display large persistent currents, of order Iq. 

In this work, Tlr[, and hence X, were calculated numerically by solving the Lippmann- 
Schwinger equation for the wire using a straight-forward generalization to three dimensions of the 
method of Nonoyama et al?^ 

The persistent current / is periodic in the flux threading the ring, fluctuates as the number of 
electrons in the ring is varied, and depends on the precise configuration of the defects. AH of this 
must be considered when specifying what is meant by the "typical" persistent current 7^^. The per- 
sistent current may be written as a Fourier series in the magnetic flux: 

/=L~l /„sin(27Un4)) 

In the present work, I^, the first Fourier coefficient, was calculated at zero Kelvin as a function of 
the number of electrons present in the ring. Its root mean square value /^^^ =^<I-^> was evaluated, 
the average <...> taken over all electron populations ranging from 1/2 to 1 electron per site.^^ /^^^ 
was then averaged over configurations of grain boundaries. This configuration average will be re- 
ferred to as Ifyp, since it is a typical value of the persistent current in an isolated ring. Since X also 
fluctuates with the Fermi energy and defect configuration, X^^, the average of X over grain bound- 
ary configurations and Fermi energy was calculated. For consistency with the above definition of 
Ifyp, the range of Fermi energies over which X was averaged was [-1.73^0, 0]; this corresponds to 
electron populations ranging from 1/2 to 1 electron per site of Hamiltonian H, for a bulk system. 



3. Results for Separable Hamiltonians: Simple Grain Boundaries 

Let us consider, to start with, rings with only "radial" grain boundaries (see left inset, Fig.l). 
In the Hamiltonian //, this means that for certain values of j (those at grain boundaries) tjjli = tgi, 
for all k and /, while all other tjj^i are equal to tQ?"^ Thus H is separable. [Similar results were ob- 
tained for non-separable Hamiltonians (for rings with small numbers of dilute grain boundaries 
running in arbitrary directions with steps on the atomic scale, and for rings with small numbers of 
atomically-rough grain boundaries) and will be discussed in Section 4.] The rings simulated had 
square wire cross-sections of NxN sites (see left inset. Fig. 1), and aspect ratios M = Z/N, for a ring 
circumference of Z sites. 

In typical thin gold films the grain size is of roughly equal to the thickness of the film ' 
which would imply that the number of grain boundaries in the ring is roughly equal to M, the aspect 
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ratio. However the grain size can be larger for some substrates and/or after heat treatment. The 
number of grain boundaries in the gold rings of Chandrasekhar et al^ ^ was not known accurately, 
but may have been significantly smaller than M. But, in the absence of more detailed informa- 
tion, it seems reasonable to examine closely models where the number of grain boundaries equals 
M, choosing the grain boundary strength parameter t^ytQ in such a way as to make the transmission 
probability X^y of the ring approximate the experimental value, X^^ ~ 0.01. This choice means that 
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in a ring with more grain boundaries, the reflection probability at each individual grain boundary 
is set to be weaker than in a ring with fewer grain boundaries so that the total transmission proba- 
bility of an electron around the ring is approximately the same. 

This is the case considered in Fig.l where the symbols o, □, and ▲ show Ify^JJ^ and 
as a function of the wire cross-section N, computed for three choices of the aspect ratio M, for radial 
grain boundaries. In each case the number of grain boundaries in the ring is equal to M. The con- 
figuration averages were computed assuming grain boundaries randomly distributed around the 
ring. In each case, the value of tgi, was chosen so that -0.01 {tgiJtQ = 0.066, 0.210 and 0.485 for 
M =1,4 and 16, respectively); see the lower right inset of Fig.l for the precise values of com- 
puted for these values of the grain boundary parameters. The straight lines in Fig. 1 are guides to 
the eye. 

Notice that for N~30, Ifyp ~Iq even though ~ 0.01. This demonstrates that a model with 
non-interacting electrons can exhibit a large persistent current ~ Iq despite the presence of scatter- 
ing that results in a short transport mean free path / ~ 0.0 IL and correspondingly small conductance 
for the ring. Note also that in Fig.l, for fixed M, Ityp/lQ increases linearly with N (there are small 
fluctuations about the straight line behavior), while slowly decreases. This behavior is in 
marked contrast with the case of random impurity scattering where Ifyp is predicted^ to be two or- 
ders of magnitude smaller (Ityp/lQ ~ X) for the same value of the conductance, and/fy^/g is predicted 
to be independent of A^. 

In Fig.l, Ifypllo also increases when the aspect ratio M (and number of grain boundaries) in- 
creases at fixed X^^, but the change of I^yp from M =4 to M =16 is small. The fact that ItypllQ in- 
creases somewhat while the number of grain boundaries is increasing is not unreasonable since the 
scattering strength of the individual grain boundaries is decreasing so as to keep the transmission 
probability X^^, of the whole ring approximately fixed. 

Since T], the number of transverse modes of the ring, is proportional to A^^, the linear depen- 
dence of Itypllo on A'^ in Fig.l resembles the result Ifyp~lQ^r\ predicted^ for ballistic (defect-free) 
rings. However here Ifyp « Iq Vt]. Also, unlike the ballistic case, here one must vary not just A'^ 
(which controls Tj) but also the circumference Z (keeping M =Z/A'^ constant) to see linear behavior. 
This is shown in Fig. 2, where IfyplI^ is plotted for rings with 4 radial grain boundaries (tgtJtQ = 
0.2 1 0) , for M =4 and 1 6, and also for fixed Z = 1 2. Increasing N (and r] ) at fixed Z results in slower 
growth of Ityp/Io than the linear behavior seen for fixed M = 4 and 16. Thus the channel counting 
arguments that yield the Vr] scaling of IfyjJlQ in ballistic rings are not adequate for rings with grain 
boundaries. For the latter, the aspect ratio M is also significant; the behavior of Ityp/Io at very small 
values of M (where the concentration of grain boundaries in the ring is relatively high) is different 
from what it is at high M where the grain boundaries are more dilute. Notice also that IjypllQ chang- 
es little when M changes from 4 to 16 in Fig.2, or when the number of grain boundaries changes 
from 4 to 16 at fixed X^^y for M =16 (compare Fig. 1 and Fig.2). 

It is of interest also to consider separately the dependence of the persistent current on the as- 
pect ratio of the ring and on the number of grain boundaries. This is addressed in Fig.3 and Fig.4. 
Fig. 3 shows examples of the dependence of IfyfJIo on the aspect ratio M for a fixed wire cross-sec- 
tion N=20 and for fixed numbers of grain boundaries. The results for 4 grain boundaries with tgi/ 
tQ = 0.210 and for 16 grain boundaries with tgfJtQ = 0.485 are shown. (These grain boundary 
"strengths" are the same as for the M=4 and M=16 results shown in Fig. 1, respectively). For small 
values of M the dependence of Ityp/Io the aspect ratio M is non-monotonic and qualitatively dif- 
ferent in the two cases shown. But at larger values of M, Ity.p/Io becomes almost independent of M, 
in agreement with the results shown in Fig.2. The fact that the behavior seen for very low values 
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of M in Fig. 3 is anomalous should not be surprising since there the grain boundaries themselves 
constitute a large fraction of the material of the ring and grossly alter its electronic structure. 

In Fig.4, •, A and ■ show Ifyp/lQ as function of the number of grain boundaries for a fixed 
aspect ratio M=4 and wire cross-section N=20, for tgiJtQ = 0.800, 0.485 and 0.210 respectively. No- 
tice that here in each case the strength the individual grain boundaries is not changed as the number 
of grain boundaries in the ring increases (in contrast to Fig.l). In Fig.4, for small numbers of grain 
boundaries, the persistent current begins to decrease exponentially with increasing number of grain 
boundaries. As the number of grain boundaries continues to increase and the grain boundaries 
cease to be dilute, ItyjJJQ begins to deviate from the exponential behavior. This deviation is stron- 
gest for the more strongly scattering grain boundaries. For weak scattering the exponential behav- 
ior persists to significantly higher numbers of grain boundaries in the ring. That the breakdown of 
the exponential behavior is principally an effect of the increasing concentration of grain boundaries 
in the ring can be seen by considering the results denoted by <> in Fig.4. These again show IfypllQ 
for A^=20 and tgiJtQ = 0.210 (as for the results denoted ■), but in this case M is not fixed but set 
equal to the number of grain boundaries, so that the concentration and the strength the of grain 
boundaries in the ring are both kept fixed when the number of grain boundaries increases. I.e. here 
the number of grain boundaries is proportional to the circumference of the ring. In this case it can 
be seen that dependence of I^p/lQ on the number of grain boundaries is closer to exponential, when 
the number of boundaries (and the ring circumference) are large. 

If the number of grain boundaries in the ring exceeds the aspect ratio M, the linear depen- 
dence of Ityp/lQ on N at constant M that is seen in Fig. 1 breaks down for small A'^. This is illustrated 
in the inset of Fig.2 where Ityp/lQ is plotted against A'^ for 4 radial grain boundaries with tgfJtQ = 
0.210 and a fixed aspect ratio M=l. The breakdown of the linear behavior at low A/^ occurs because 
if the aspect ratio and the number of grain boundaries are both kept fixed, the grain boundaries 
themselves take up an increasingly large fraction of the material of the ring as A'^ decreases. As has 
already been seen for low aspect ratios M in Fig. 3, this changes the underlying electronic structure 
of the ring very strongly and modifies the behavior of the persistent current. Clearly the gold rings 
of Chandrasekhar et al}^ where the grain boundaries most likely comprise less than 0.1% of the 
material of the ring, are in a dilute regime where this effect does not occur. But it becomes a major 
concern when modelling rings with non-separable Hamiltonians, since only quite small systems of 
that type can be simulated numerically, which severely limits the number of grain boundaries in 
rings that can be modelled if the grain boundary concentration is to be kept low. 

4. Results for Non-Separable Hamiltonians 

The above results are for Hamiltonians with radial grain boundaries. It is necessary to estab- 
lish whether they also apply to systems with grain boundaries running in arbitrary directions. The 
Hamiltonians of such systems are non-separable, and their numerical study is much more difficult, 
requiring the use of Lanczos algorithms to calculate the electronic energy spectra. Some repre- 
sentative results are given in Fig. 5. Here, for the results denoted • and o, the number of grain 
boundaries is equal to the aspect ratio M, as in Fig. 1, and t^ytQ is chosen so that -0.01. The 
grain boundaries considered are planes slicing all the way through the ring, as is illustrated in the 
right inset of Fig.5. Planar grain boundaries are realistic, being favored by free-energy consider- 
ations during sample preparation. But in the simulations the "planes" were in fact stepped interfac- 
es on the atomic scale of the lattice Hamiltonian H. Configuration averaging was performed over 
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the possible orientations of such grain boundaries, but excluding boundaries parallel to the (3 or y 
(or both the |3 and y) hopping directions in the Hamiltonian. This exclusion ensured that averaging 
was only over non-separable Hamiltonians. These results in Fig.5 are similar to those in Fig. 1 : 
slowly decreases with increasing A^. Ifyp/Io increases linearly with for one grain boundary with M 
=1 (•). For two grain boundaries with M=2 ( o), the behavior is also consistent with linear growth 
of Ifyp/lQ with A^, but the fluctuations are larger, and the accessible values of A'^ smaller. (The fluc- 
tuations are a finite size effect; they do not disappear when more grain boundary configurations are 
included in the average Ifyp.) These results indicate that the behavior described above for radial 
grain boundaries is not specific to separable Hamiltonians. The most obvious difference between 
the above results for the systems with radial and non-radial grain boundaries is that in the latter 
case Ifyp is smaller by a factor ~ 2, for the same value of X^^. 

A feature of the results for 2 grain boundaries in the M=2 structure that are denoted o in 
Fig.5, is that at low A'^ where the concentration of grain boundaries in the crystal is highest, Ifyp/Io 
is almost independent of N (for N <6). A similar but less pronounced tendency for IfyjJiQ to level 
out at low A'^ can also be discerned in the results • for the M=l structure with one grain boundary. 
If the grain boundaries in the ring are made more concentrated, this tendency of IfypllQ to level out 
at low N becomes stronger. This is illustrated by the results denoted ♦ in Fig.5 which show IfypllQ 
for 2 grain boundaries in an aspect ratio M =\ structure, with tgi,ltQ = 0.076. Here Ifyp/Io is almost 
independent of A'^for A^<15, and then shows an increase, but it is not clear whether the linear growth 
of Ifyp/lQ with A'^ is recovered at large N or not, because of the fluctuations that IfypllQ shows as a 
function of N and because of the limited range of A'^ that is accessible numerically. 

Notice that the breakdown of the linearity of ItyjJlQ with A^ over the wide range of A'^ seen for 
the 2 grain boundary M=l system (♦ ) does not occur in the 2 grain boundary M=2 system ( o) or 
in the 1 grain boundary M=\ system (•). Thus neither the number of grain boundaries alone nor 
the value of M alone controls this effect. An important factor appears to be the grain boundary con- 
centration. These results (like those presented in Fig.2, Fig.3 and Fig.4, and the discussion at the 
end of Section 3) demonstrate the importance of simulating systems in which the grain boundaries 
are dilute, as they are in the real gold rings. However, for models with larger values of the aspect 
ratio M, the numerically accessible range of A'^ shrinks very rapidly. This, together with the fluctu- 
ations that IfypllQ shows as a function of A^, makes it very difficult to establish how IfypllQ scales 
with A'^ for large numbers of dilute non-radial grain boundaries. 

Although the grain boundaries discussed above correspond to non-separable Hamiltonians, 
are stepped on the atomic scale and are randomly located and oriented in the ring, each of them has 
a perfectly regular microscopic structure (except where two grain boundaries cross). Recent high 
resolution electron microscopy studies^^ have yielded very detailed and clear images of grain 
boundaries in gold that are perfectly ordered on the atomic scale. However, earlier work indi- 
cated some grain boundaries in gold to be very well ordered, but others to exhibit some atomic 
scale disorder in a region along the grain boundary a few atoms thick. Since the micro structure of 
the grain boundaries of the gold rings of Chandrasekhar et al. ^ ^ was not measured, it is unclear 
which type of grain boundary was more typical, and it is of interest to simulate atomically rough 
grain boundaries as well. 

The results of simulations for a model of a ring with an atomically-rough grain boundary, are 
shown by the + symbols in Fig. 1 for M = 2. This grain boundary contains a central layer identical 
to one of the radial grain boundaries described in Section 3, but surrounded by two disordered lay- 
ers. The disordered layers are constructed by making 50% of the hopping matrix elements in the a 
direction (at randomly chosen sites) in the lattice layers adjacent to that of the central layer also 
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weak and equal to tgi^. The results for this Hamiltonian are again similar to those for the simple 
radial grain boundaries shown in Fig. 1 . The effect on the (configuration averaged r.m.s.) persistent 
current of roughening the grain boundary in this way is similar to that of tilting the grain boundary, 
i.e., making it non-radial. 

5. Comparison of Grain Boundary and Random Defect Models 

Why do these grain-boundary models yield the combination of large persistent currents and 
small conductances that is observed experimentally, while random defect models do not? This is 
addressed in Fig.6. Here Ifyp/lQ and are results for M =1 (one radial grain boundary) taken from 
Fig. 1. and X.^ are the persistent current and transmission defined in the same way as 7^,^ and K^y, 
but with the same "defects" (terms in H for which tjli = tgi^) scattered randomly through the system 
instead of forming a continuous grain boundary. The striking feature is that while is smaller 
than Xf. by a factor -40, I^yp is smaller than 7^ by only a factor ~2. That is, collecting the random 
defects into a coherent grain boundary forming a barrier extending all the way across the path of 
the current, reduced the transmission A, (conductance of the ring) far more than it did the persistent 
current. 

This clarifies why large persistent currents occur for low conductances in the grain boundary 
model, but not in random defect models. The physics underlying this numerical example is that the 
mean free path is a transport property, while the persistent current is an equilibrium effect. The 
persistent current is given by the flux-derivative of the total electronic energy of the ring. The total 
energy (and hence the persistent current) is not very sensitive to the arrangement of the defects, and 
not disturbed much if defects form into a grain boundary, i.e., a continuous "barrier" across the 
ring. But transport (the conductance or X) is well known to be greatly affected by such a barrier. 
Analogous (and more extreme) differences between the behavior of transport coefficients and 
equilibrium quantities are commonly associated with percolation phenomena where the topology 
of the sample is crucial. 

This physical argument suggests that the present results should not be very sensitive to the 
shape of the grain boundary, so long as it forms a continuous barrier across the ring. For example, 
as has been shown above for a ring with one grain boundary, making the grain boundary rough on 
the atomic scale does not change the predictions of the model qualitatively. 

It should be stressed that it is the conductance rather than the persistent current that is sensi- 
tive to the arrangement of defects in Fig.6. Thus the present work suggests that the explanation of 
the experiments of Chandrasekhar et al. ^ ^ has at least as much to do with the behavior of the con- 
ductance of the ring as with that of the persistent current, and that the conductance aspect of the 
problem deserves much more attention than it has received in previous theoretical studies. 

It is also clear from the above discussion of Fig.6 that, in order to achieve a low electron 
transmission probability (X ~ 0.01), the random impurity model requires a much larger number of 
atomic defects to be present in the ring than is required in the grain boundary model (if the indi- 
vidual atomic defects are similar). This in turn depresses the value of the persistent current that the 
random impurity model predicts for a given conductance. Since the resistance of gold films is 
mainly due to scattering by grain boundaries, ' the high concentration of defects that is required 
by random impurity models in order to obtain low conductances is a deficiency of those models. 
That is, the random defect models assume the gold rings to be much "dirtier" than they really are. 
For example, computer simulations of persistent currents in rings with random impurities have tra- 
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ditionally employed Anderson models in which every site is an impurity site and the spread of ran- 
dom site energies W is not very much smaller than the tight binding energy band parameter ^q- 

6. Comparison with Experiment 

The amplitudes of the persistent currents measured by Chandrasekhar etal^^ ranged in mag- 
nitude from O.S/q to 2.0Iq over the three rings for which results were reported. For these rings, N ~ 
300, so that the largest wire cross-sections N considered in Fig.l and Fig. 5 are about an order of 
magnitude smaller than in the experiments. The aspect ratios of the experimental rings were in the 
range M ~ 100 to 200. The numbers of grain boundaries in the rings were unknown, but the trans- 
mission probabilities of the rings could be estimated from the known conductances of similar sam- 
ples to yield X ~ 0.01, as discussed in Section 2. If one assumes that the rings contain any given 
number of grain boundaries, one can choose the grain boundary model parameter tgjJtQ in such a 
way that the electron transmission probability of the model ring matches that in the experiments, 
and this was done in the calculations presented in Fig.l and Fig.5. 

One can extrapolate the results of the simulations reported in Section 3 for rings with radial 
grain boundaries to systems of with the experimental values of N and M, by using the linearity of 
Ityp/Io with N demonstrated in Fig. 1 and the fact that for fixed N and fixed numbers of grain bound- 
aries, Ifyp/lQ is independent of M at large M, as shown in Fig.3. This yields /^^ ~ 10 /q for rings with 
the experimental dimensions, if the rings contain 4 or 16 grain boundaries and the ring conductance 
matches experiment, i.e. X ~ 0.01. It should be noted that in the models with radial grain bound- 
aries, for large numbers of grain boundaries and large M, I^ypllQ is insensitive to both M and the 
number of grain boundaries ifX is held fixed. Thus the prediction Ip^p ~ 10 /q for the radial grain 
boundary models is insensitive to the number of grain boundaries that one assumes the ring to con- 
tain, if the ring's N, M and conductance are required to match their experimental values. 

For non-radial grain boundaries, one can extrapolate the linear behavior of the results seen 
in Fig.5 for M = 1 structures with one grain boundary, and for M=2 structures with two grain bound- 
aries, to the experimental value of A^, while holding M fixed. This yields I^yp ~ 3Iq if A'^ and the con- 
ductance of the ring take the experimental values. For 2 grain boundaries, the results of Fig.5 show 
that Ifyp IIq increases when M increases from 1 to 2. Thus it is not unreasonable to expect rings with 
one or two non-radial grain boundaries, and having the experimental values of A^, M and the con- 
ductance, to exhibit large values of the persistent current Ifyp > 3Iq, that are quite similar to those 
observed by Chandrasekhar et al.^^ However, as discussed in Section 4, for large numbers of non- 
radial grain boundaries the important dilute regime (of laxge N and M) is not readily accessible to 
numerical investigation; a definitive treatment of that regime is beyond the scope of the present 
work. 

It should also be emphasized that both the separable and non-separable Hamiltonians con- 
sidered in this article represent very simple models of grain boundaries. Complete ab initio calcu- 
lations of more realistic grain boundary structures, of the associated electronic eigenstates and of 
the persistent currents and transport in rings containing them would clearly be of interest. 

To summarize: There is good experimental evidence that grain boundary scattering is the 
dominant electronic scattering mechamsm in thin gold films. ' The present simulations have 
demonstrated that exact solutions of some simple models that include grain boundary scattering 
exhibit the combination of large persistent currents and small conductances that is observed exper- 
imentally in gold rings, ^ ^ but cannot be explained by random impurity models with non-interacting 
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electrons. The random impurity models yield persistent currents that are too small by factors of 30 
to 150.^^ However, it is unclear how many grain boundaries the gold rings of Chandrasekhar et 
al. ^ ^ actually contained, and whether more realistic models or models with many non-radial grain 
boundaries can account for the experimental data as well as the simple models discussed above do, 
is yet to be determined. 

7. Proposed Experimental Tests 

It is clear that direct experimental tests of the predictions of the grain boundary scattering 
explanation of the persistent current measurements of Chandrasekhar et al. ^ ^ are needed, and some 
are suggested below: 

1) The experimental data reported by Chandrasekhar et al}^ exhibited a sample-to- sample 
variation of the size of the observed persistent current (///q) by a factor of ~ 7, over the three sam- 
ples. This is consistent with the idea that scattering by grain boundaries (whose number varies ran- 
domly from sample to sample) controls the magnitude of the persistent current; see, for example, 
Fig.4. However, similar variations in the magnitude of the persistent current are also found to occur 
in simulations if the number of grain boundaries is kept fixed but their configuration changes. (Re- 
call that the numerical values of I^yp presented in this paper are grain boundary configuration aver- 
ages of r.m.s. values, as defined in Section 2.) It would be of interest to establish experimentally 
whether there is a correlation between the magnitude of the persistent current and the number of 
grain boundaries that the sample contains, by performing persistent current measurements fol- 
lowed by spatially resolved electron diffraction measurements on a series of samples. But in such 
a study it would be important to bear in mind the sensitivity of the persistent current to the details 
of the grain boundary configuration; comparisons between a small number of samples with similar 
numbers of grain boundaries would not be significant. It would also be very desirable to measure 
the conductances of the same samples directly instead of relying on conductance measurements on 
similar samples, as has been done in the past. 

2) It should be possible to fabricate single crystal gold rings that do not contain any grain 
boundaries at all. These should exhibit significantly larger persistent currents than those that have 
been observed to date. 

3) If single crystal rings of gold alloys can be fabricated, these should exhibit weak persistent 
currents for low conductances, since the random impurity theories^' should apply to such sys- 
tems. 

1 '~) 

4) Semiconductor rings such as those of Mailly et al. do not contain any crystal grain 
boundaries in the conducting region. The quasi-ballistic rings of this type that have been studied to 
date have shown no evidence of any anomalous behavior; the magnitudes of the persistent currents 
observed in them are consistent with the predictions of standard, non-interacting electron theories, 
which for this case predict that /^^ ~ Iq since the number of channels is small and / ^ L. By intro- 
ducing impurities directly into the region occupied by the conducting electrons it should be possi- 
ble to fabricate diffusive semiconductor rings for which l«L but which are still free of grain 
boundaries. It would be very interesting to establish whether or not large persistent currents ~ Iq 
occur in such samples. However, it is necessary to measure both the persistent currents and con- 
ductances of such samples for this test to be meaningful. It should also be stressed that a truly dif- 
fusive system, preferably with point defect type scatterers, is required for this test to be convincing 
~ Using charged random impurities with long-range interactions can result in a few electronic 
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modes following a preferred path through the ring and thus yield misleading results. 

5) Although the semiconductor rings are free of grain boundaries, it is possible to impose an 
electrostatic potential barrier (the analog of a grain boundary) across such a ring by means of a gate. 
In such a system it should be possible to measure both the transmission probability of electrons 
through the barrier and the persistent current in the ring. It is predicted^^ that in the regime of quan- 
tum tunneling of the Fermi electrons just below the top of the barrier, large persistent currents ~ Iq 
can coexist with small electron transmission probabilities « 1. It would be of interest to observe 
this semiconductor analog of the large persistent currents that coexist with low conductances in the 
small gold rings. 

I wish to thank C. Barnes, B. Heinrich, B. L. Johnson, D. Loss and R. A. Webb for interesting 
discussions. This work was supported by the Natural Sciences and Engineering Research Council 
of Canada. 
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Figure Captions 

Fig.l Normalized persistent current IfypllQ and average transmission vs. N, for radial grain 
boundaries, are shown as o, □ , and ▲ for ring aspect ratios M =1,4 and 16 respectively. Number 
of grain boundaries = M. Straight lines are guides to the eye. + symbols show results for an atom- 
ically rough grain boundary, for M=2 and tgjJtQ = 0.047. Left inset: Schematic of ring with 4 radial 
grain boundaries (shaded). 

Fig.2 IfypllQ vs. A'^ for 4 radial grain boundaries. □ , +, ♦ are results for fixed aspect ratios M = A 
and 16, and a fixed circumference Z =12, respectively. Inset: IfypHo vs. for 4 radial grain bound- 
aries with tgiJtQ = 0.210 and a fixed aspect ratio M = 1. Note the breakdown of linearity of Ifyp/lQ 
with N at low A^^ where the grain boundaries themselves constitute a significant fraction of the ma- 
terial of the ring. 

Fig.3 Ityp/Io vs. aspect ratio M for a fixed wire cross-section N = 20 and for fixed numbers 4 and 
16 of radial grain boundaries. 

Fig.4 •, ▲ and ■ show Ifyp/lQ vs. the number of radial grain boundaries for a fixed aspect ratio 
M=4 and wire cross-section N-20, for tgiJtQ = 0.800, 0.485 and 0.210 respectively. O also show 
Ifyp/Io foi" ^=20 and tgiJtQ = 0.210 (as for the results denoted ■), but with the aspect ratio M set 
equal to the number of grain boundaries, so that the concentration of grain boundaries in the ring 
is kept fixed. The straight lines are guides to the eye. 

Fig.5 IfypllQ and vs. A^, for non-radial grain boundaries. • and o are results with the number 
of grain boundaries equal to the aspect ration M, for aspect ratios M =1 and 2 and tg^tQ = 0.060 and 
0.075 respectively. The results denoted ♦ show Ityp/lQ for 2 grain boundaries in an aspect ratio M 
=1 structure, with tgiJtQ = 0.076. Error bars are the statistical uncertainty due to averaging over fi- 
nite numbers of grain boundary configurations. Straight lines are guides to the eye. Right inset: 
Schematic of 2 non-radial grain boundaries. 

Fig.6 Comparison of random defect and grain boundary models. See text. 



